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ABSTRACT 

On 13 February 2009, a coronal wave - CME - dimming event was observed in quadrature by the 
STEREO spacecraft. We analyze this event using a three-dimensional, global magnetohydrodynamic 
(MHD) model for the solar corona. The numerical simulation is driven and constrained by the 
observations, and indicates where magnetic reconnection occurs between the expanding CME core 
and surrounding environment. We focus primarily on the lower corona, extending out to 3i?©; this 
range allows simultaneous comparison with both EUVI and COR1 data. Our simulation produces a 
diffuse coronal bright front remarkably similar to that observed by STEREO/EUVI at 195 A. It is 
made up of two components, and is the result of a combination of both wave and non-wave mechanisms. 

The CME becomes large-scale quite low (< 200 Mm) in the corona. It is not, however, an inherently 
large-scale event; rather, the expansion is facilitated by magnetic reconnection between the expanding 
CME core and the surrounding magnetic environment. In support of this, we also find numerous 
secondary dimmings, many far from the initial CME source region. Relating such dimmings to 
reconnecting field lines within the simulation provides further evidence that CME expansion leads 
to the "opening" of coronal field lines on a global scale. Throughout the CME expansion, the coronal 
wave maps directly to the CME footprint. 

Our results suggest that the ongoing debate over the "true" nature of diffuse coronal waves may 
be mischaracterized. It appears that both wave and non-wave models are required to explain the 
observations and understand the complex nature of these events. 

Subject headings: MHD - Sun: corona - Sun: Coronal Mass Ejections (CMEs) 
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1. INTRODUCTION 

The physical nature of EUV coronal wav es has been 
ambiguous ever since their discovery in 1996 (fpere et al.l 
ll997HMoses et al.lll997llThompson et al.lll998D. with the 
Extre me Imaging Telescope (EIT: iDelabo udinicr e et al.l 
119951 ) o n board the Solar an d Hcliospheric Observatory 
(SOHO: [Domingo et aLlll995l ). Following the launch of 
the Solar Terrestr ial Relations Observatory (STEREO; 
iKaiser et~a l. 2008), coronal waves are now observed with 
increased temporal and spatial resolution, and have now 
been seen from two persp ectives using the E xtreme 
Ultra- Violet I magers (EUVI: [Wuelser et al.l l2004l). Using 



SOHO data, iBiesecker et al. ~ (f2002l ) showed (after cor 



recting for observational biases), that every EIT wave in 
their study could be credibly associated with a coronal 
mass ejection (CME). 

In the past decade, there has been ongoing debate 
about the physical nature of the coronal wave bright 
front. Many models have been developed and they can be 
split into two main groups. The first group of models de- 
scribes coronal waves as a pure MHD wave, either a: fast- 
nrr; : — n i i rrTi ' : — n i »»„» i 



mode 



1999; Klasscn et~aT1 1200C 



Vrsnak et all 120021 iCliver et all l20Q4t iWarmuth et alj 


20049), slow-mode ( 


Krasnoselskikh & Podladchikova 


2007: Wans et al.l 12009 


Y, or a slow-mode solitary wave 



In many respects, a wave 
model is appropriate. Fast-mode solutions are consis- 
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tent with propagation across the magnetic field; slow- 
modes can account for the large intensity enhancements 
and some of the lower measured velocities; and solitary 
waves are consistent with single-pulse structure and co- 
herence over large distances. In these models, the wave 
front may be freely propagating (generated by an ini- 
tial pressure pulse), or piston-driven where the wave is 
cons tantly supplied by ener gy from the expanding CME 
(e.g. lVrsnak & CliverlhoOST ). 

The second group of models describe the bright front 
as being related to the actu al expansion of the CME, in a 
non-wave contex t . Work by[D clannc e fc Aulanierl (|1999|) ; 
iDelanned (20001 ) ; lDelannee et all (|2008D attribute the re- 
gions of enhanced emission at low altitude to compression 
of the plasm a between stable flux domains. iChen et alJ 
(|2002t 120051 ): IChen! (|2009l ) suggest that the expansion 
of the magnetic field during the CME lift-off gradually 
"opens" the overlying magnetic fi eld, compressing the 
plasma in the legs of the CME. And I Attrill et all (|2007af) 
propose magnetic reconnections between the expanding 
CME and favorably oriented surrounding quiet Sun mag- 
netic field as an explanation for the patchy bright front. 
Although the connection with CMEs factors strongly 
in many of these models, becaus e the exact relation- 
ship between diffuse coronal waves (|Biesecker et al.l l2002: 
!Vrsnakl2005D and CMEs is not well understood, the phys- 
ical nature of the bright front has remained inconclusive. 

Space-based white-light coronagraph observations of 
the low corona are available from STEREO/COR1. 
However, the field of view only starts at 1.3 Rq, and 
the domain where coronal waves are observed, below 
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200 Mm (e.g. iWarmuth et all 


2004al 120051: IVrsnak et all 


1200511 Wills-Davev et al.ll2007l 


Patsourakos et alj 2009f) is 



blocked by the coronagraph occulting disk. Thus, a com- 
bination of EUVI and COR1 observations, coupled with 
numerical simulations, is required to develop a coherent 
picture of the early stages of CME evolution in the low 
corona. 



In this paper, we present an analysis of observations 
combined with a global numerical simulation of a coro- 
nal wave event observed by STEREO on 13 February 
2009. This is the first numerical simulation of an EUV 
coronal wave based on this observation in quadrature. 
We use a global MHD model for the solar corona that 
is driven by real magnetogram data, and we drive the 
CME in a realistic way that matches the observations. 
Numerical studies of coronal waves have been carried 
out previously. H owever, these model s were either two- 
dimensional (e.g. IPomoell et aT1l2008D . considered only 
the expanding flux tube, neglecting interaction with the 
surroundings (jDelannee et al.1 [2008). simulated the in- 
terac tion of coronal waves with only a local active re- 
gion (|Uchidal[l97l IQfman fc Thompson! [20021) . or they 
drove the coron al wave by a pressure pulse (|Wanell2000t 
IWu et al.ll200ll) . 

By constraining the simulation as much as possible 
with the observations, we obtain a full picture of the 
three-dimensional evolution of the coronal magnetic field 
during the eruption, including those regions that are not 
observed by white- light coronagraphs. In this way, we 
hope to shed some light on the nature of coronal waves, 
their relationship to CMEs, and their theoretical descrip- 
tion. 

We describe the observations of the CME event in Sec- 
tion [2 and the numerical simulation in Section [3] We 
present the results in Section [4] and discuss the impli- 
cations for the various descriptions of coronal waves in 
Section \5\ We conclude our findings in Section [5] 

2. OBSERVATIONS 

The coronal wave - CME - dimmings event on 13 
February 2009 occurred when the STEREO spacecraft 
had a separation angle of 91°. The Sun was in the very 
quietest part of its cycle, just at the start of the rise 
phase of solar cycle 24. NOAA active region 1012 was 
the only active region on the solar disk. Thus the major- 
ity of the surrounding environment was quiet Sun, with 
a low-latitude coronal hole to the East of the active re- 
gion (see top panels Figure [T]). Pre-eruption, the active 
region hosted a forward "S" sigmoid. The Sun produced 
a CME that was seen by spacecraft A right on the East 
limb, and by spacecraft B expanding from the center of 
the disk. 

2.1. STEREO/COR1 

The COR1 instruments are internally occulted coron- 
agraphs and obs erve the inner solar coro na in white light 
from 1.3-4 R ([Thompson et al.1l2003h . Base difference 
images (where a pre-event image at 05:45 UT, is sub- 
tracted from all subsequent images) of COR1-A data, 
are shown in the left panels of Figure The COR1 im- 
ages were taken with a temporal cadence of 10 minutes, 
and have a pixel size of 7.5". CORTB also observed the 
CME as a halo event, first becoming apparent beyond the 



occulting disk at 06:55 UT. The non-differenced COR1-A 
data (not shown) show a helmet streamer located north 
of AR 1012, and open streamers to the south. 

2.2. STEREO/EUVI 

The EUVI imagers observe the Sun out to 1.7 R Q , 
and produce 2048x2048 images with a pixel size of 1.6" 
([Wuelser et al.ll2004h . We analyze the 195 A EUVI im- 
ages, which have a temporal cadence of 10 minutes. 
Base difference images from EUVI-B are shown in Fig- 
ure [3J The base difference images are produced by first 
compensating for the s olar rota tion usin g SolarSoft's 
drot_map routine (http://www.lmsal.com/solarsoft), so 
that all images are de-rotated to the pre-event image time 
at 05:15 UT. Then, the pre-event image is subtracted 
from all subsequent images. Base difference images high- 
light real intensity changes, with bright areas showing an 
increase in emission with respect to the pre-event data. 

2.2.1. Coronal dimmings 

As well as showing the coronal wave bright front, 
base difference images also show depletions in intensity, 
known as "coronal dimmings" . These are regions where 
the plasma density has dramatically decreased due to 
plasma evacuation along the "opened" magnetic field, 
usual l y occurring during an eruption (e.g. iHudson et al.l 
[19961: |H arm fc Sterling! I200II : IHarra et al.N2007l) . Vari- 
ous works have shown that the core coronal dimming 
regions (i.e. black regions, bottom panels of Figured]), 
located on either side of the bright post-eruptive arcade 
(PEA) mark the footpoints of the erupting flux rope (e.g. 
i Webb et al.l2000tlMandrini et alj2005HCrooker fc Webbl 
l2006t lAttrill et alj|2006t iMcIntosh et al J 120071) . The in- 
tensity drop in core cor onal dimmings is substantia l (typ- 
ically ~ 40 - 60%; e.g. lChertok fc Grechnevll2005l ). 

We process the EUVI base difference data using the 
automatic dimmings detect ion algorithm described in 
lAttrill fc Wills-Davevi (|2009D . Figure U shows the out- 
put from this algorithm. In addition to the core dim- 
mings near to the PEA, secondary dimmings are also de- 
tected which develop remote from the active region and 
are spread across the solar disk. These secondary dim- 
mings are more subtle than the deep, core dimmings and 
are not easy to identify by eye in the base difference data. 
We discuss the secondary dimmings further in t|5.4l 

2.2.2. Persistent brightenings 

The base difference images also show relatively concen- 
trated, persistent brightenings at the edge of the deep, 
core dimmings. Examination of the location of coronal 
holes with respect to these brightenings shows that two 
brightenings (marked "A" and "B" in Figure [1} are sit- 
uated along the boundary of a low-latitude coronal hole 
extending to the East of the active region. Additionally, 
brightenings "B" and "C" are located at the edge of the 
core dimmings. 

3. NUMERICAL MODEL 

In order to simulate the CME event, we use the Solar 
Coro na (SC) model develop e d at the University of Mich i- 
gan (|Roussev et al.l l2003bl : iCohen et all [20071 l2008bl) . 
The model is based on the global MHD BATSRUS code 
(|Powell et al. 1999) and is part of the Space Weather 
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Modeling Framework (SWMF) (jToth et al.ll2005f) . The 
model solves the set of MHD equations on a non-uniform 
Cartesian grid and is designed with highly efficient par- 
allel architecture. Here, we describe the model briefly. 
We refer the reader to the references above for a more 
detailed description. 

The ambient solar wind conditions in the model are ob- 
tained under the assumption that the source of energy re- 
quired to power the solar wind is the change in the poly- 
tropic index, 7 in a non-polytropic medium. The numer- 
ical procedure is as follows. First, the potential magnetic 
field is calculated using high resolution SOHO/M ichelson 
Doppler Imager (MDI; Scher rer et alj I1995D magne- 
togram synoptic maps (http://soi.stanford.edui. Second, 
this potential field distribution is used to calculate the 
distribution of the terminal solar wind speed, u wsa as a 
function of the flux tube expansion fac tor, f s , based on 
the Wang-Sheeley-A r ge (W SA) model (jWang fe Sheelevl 
fl990t lArge fc Pizzol l200l . Third, the photospheric 
boundary conditions for 7, and the terminal speed, u wsa , 
are related by tracing the total energy (Bernoulli Inte- 
gral) along the flux tubes. The spatial distribution of 7 is 
then specified as a radial function of the photospheric val- 
ues, and the MHD equations are solved self-consistently 
until a steady-state with a wind solution is obtained. Fig- 
ure [5] shows the steady-state values of number density, 
n, magnetic field strength, B, temperature, T, plasma (3, 
sound speed, C s , and Alfven speed, va respectively, at 
height of r — l.li? Q . 

In order to drive the CME, we superimpose 
an unstable, sem i -circu lar flux rope based on the 
iTitov fc Demoulinl (1999) model, on top of the ambient 
solution ( Roussev et al.ll2003af ). We match the location 
and orientation of the flux rope to those of the source 
active region and its inversion line as they appear in the 
magnetogram data. The free energy of the CME is ob- 
tained by a prescribed toroidal current; we modify the 
magnitude of this current to match the observed CME 
speed. We would like to stress that even though the CME 
initiation method is not based on actual photospheric 
motions, it has been successful in mimicking the CME 
conditions once it is already propagat i ng and expanding, 
(iLugaz et al.ll2007t iCohen et al.|[2008allManchester et all 
2008). Since we are interested in the development of the 
CME after it has already been initiated, this model is 
appropriate for our study. 

We run the simulation in a Cartesian box of 2Oi? x 
20-R© x 2Oi?0 , with 9 levels of grid refinement around the 
solar surface. The grid size around the active region and 
up to a height of about 3i? is of the order of l/5Oi? . 
The ambient coronal conditions are driven by MDI mag- 
netogram data for Carrington Rotation 2080. The MHD 
simulation was performed using the Pleiades cluster at 
NASA's Advanced Supercomputing (NAS) center. 

4. RESULTS 

4.1. CME Expansion 

Prior to a detailed analysis of the magnetic field 
three-dimensional evolution, we validate the timing of 
the simulated CME by comparing its propagation with 
STEREO/COR1 observations. This validation is re- 
quired for any further assumptions about the dynamic 
evolution and its relation to the observed coronal waves. 



Fortunately, the 13 February 2009 event was observed 
by both STEREO-A and STEREO-B simultaneously, in 
quadrature. 

Figure [2] shows a side view of the CME propagation, 
comparing STEREO-A COR1 base difference white-light 
images with synthetic base difference white-light images 
generated from the simulation domain for 5:55, 6:15, and 
6:35 UT. We would like to stress that both real and simu- 
lated data sets have been processed in the same manner, 
and the scale of the images has been chosen so that it 
provides the best display. 

The comparison of the simulation (center column) with 
the base difference images (left column) is used to verify 
the structure of the CME. Comparison with the run- 
ning difference images (right column) is used to verify 
the lateral extent. We emphasise that the outer shell 
of the CME is a very subtle feature, and is only really 
evident when successive frames are switched back and 
forth. Thus the reader is encouraged to examine the 
COR1 running difference movie (CORljrdif f .mov), us- 
ing the white arrows in the right column of Figure [2] as 
a reference. One can see that the simulated CME front 
and the global structure matches well to the observed 
CME. 

4.2. Coronal Wave Expansion 

Figure [3] compares the simulation results with 
STEREO-B EUVI data for the period 5:25-6:35 UT. The 
left panel of each pair shows an EUVI base difference 
image, while the right panel shows base difference im- 
ages of the simulated mass density at a height of 1.1 i? 
(about 70 Mm) . We choose this height as it is consistent 
with calculations of th e coronal wave bright fr ont alti- 
tude from EUVI data (|Patsourakos et ai1l2009f) . as well 
as with previous height estimates from SOHO /EIT data 

We mark the leading edge of the bright front in each 
base difference EUVI-B image with a dashed white cir- 
cle. We emphasize that this circle has been drawn by eye, 
and simply indicates the maximum extent reached by the 
coronal wave in a given image. (The reader is encour- 
aged to view the supplemental movies provided with this 
paper as well: 195b_dif f .mov, densityf rontdif .mov). 
Figure [3] shows that the expansion of the simulated wave 
front is in a good agreement with the observed one. Devi- 
ations are probably due to the fact that the actual EUVI 
emissions are not simply a representation of the mass 
density, but a rather complicated, integrated function of 
it. 

The simulated bright front in Figure[3]becomes broader 
as it expands further from the active region. This is con- 
sistent with obser ved properties of coronal waves from 
SOHO /EIT data dDere et alJ[l997t iKlassen et ai] |205ot 
iPodladchikova fc Berghm ans 20Q5|). Areas of decreased 
mass density (corresponding to the core coronal dimming 
regions) can also be identified in the simulated data. We 
note that both the real and simulated bright fronts have 
an increasingly patchy, diffuse nature as the front ex- 
pands away from the active region. We measured the ex- 
pansion of the leading edge of the bright front in running 
difference data from 05:45 - 06:15 UT, which expands 
with an average velocity of f» 260 kms -1 . For compar- 
ison, the Alfven velocity from the simulation (Figure [5]) 
is w 200 kms" 1 . 
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4.2.1. Two-component bright front 

The simulations show a predominantly two-level bright 
front. Figure [3] shows that the highest intensity is con- 
centrated in patches located within a weaker, broaden- 
ing front. Although these two components expand in 
a coupled manner during the early phase of the CME, 
during the later frames, the simulated wave front is in- 
creasingly dominated by the weaker intensity component, 
which continues to expand (densityf rontdif .mov; 
densitysidedif .mov). 

The observations become increasingly noisy as the 
event progresses and the coronal wave becomes more and 
more difficult to detect. In EUVI-A base difference data, 
06:25 UT is the last frame where the bright front to the 
west of the active region is discernable (195a_dif f .mov). 
In the EUVI-A base difference simulation, the higher in- 
tensity patch to the west of the active region is visible 
until 06:35 UT (densitysidedif .mov). 

In the EUVI-B base difference data, the furthest 
reaches of the bright front (most obviously near to the 
south polar coronal hole) can be identified until 06:55 
UT. The bright front is approximately stationary at this 
time (195b_diff .mov). The EUVI-A and B base dif- 
ference simulations show that isolated higher intensity 
patches still exist at 07:15 UT. Near to the south po- 
lar coronal hole, the higher intensity patch is located at 
the same place from 06:45 UT. Near to the north polar 
coronal hole, a new higher intensity patch develops from 
07:05 UT, also remaining at the same location. 

4.3. Comparison of CME and Coronal Wave 
Expansion 

Figure [6] compares snapshots at 06:05 UT from COR1- 
A, the simulated CME, and EUVI-A, all scaled to the 
same size. The COR1-A and EUVI images are running 
differ ence images, the simulation is a base difference im- 
age. iPatsourakos &: Vourlidasl (120091) show a fit of the 
3D CME model of |Tliermsjeii_e±alJ (2006, 2009) to the 
COR1-A data for this event at 06:05 UT. This leads them 
to determine an extent for the CME in the low corona 
that is much too small to match the coronal wave in the 
corresponding EUVI data, since the 3D model is essen- 
tially made up of a spherical bubble attached to a conical 
leg. We note that a sim ilar approach was also used in 
IPatsourakos et alJ 1)20091 ). In both papers, the authors 
interpret the apparent misfit between the CME model 
extension in the low corona and the EUVI coronal wave 
as evidence that the coronal wave and CME are differ- 
ent structures and conclude that the coronal wave is a 
fast-mode MHD wave. 

Our simulation results at 06:05 UT are shown in the 
middle panel of Figure [H The simulation gives informa- 
tion about the low corona below 1.3 Rq (200 Mm), the 
region obscured by the occulting disk in the COR1 data. 
Comparison of the middle and right panels of Figure [5] 
show that the extension of the CME in the low corona 
maps very well to the coronal wave in the EUVI base 
difference data. 

In particular, the simulation results show a secondary 
cavity located to the north of the main CME cavity 
(marked by the white arrow in the middle panel of Figure 
[5]). Comparison with the corresponding EUVI base dif- 
ference data shows secondary coronal dimmings develop- 



ing at the same location (right panels, EUVI (A), Figure 
0]). Despite the lack of spectral diagnostics for secondary 
dimmings, this correlation between the secondary CME 
cavity and the secondary coronal dimmings is consistent 
with plasma evacuation. 

A time-series movie of the simulated COR1-A white- 
light data (CORljointl .mov) shows that the secondary 
cavity expands and merges with the main CME cavity, so 
that the low corona is really "opened" to a large lateral 
extent. This analysis demonstrates the important role 
of the simulation in developing an understanding of the 
true lateral extent of the CME in the low corona. 

In i j4.2.1l we noted that the higher intensity patches of 
the coronal wave front no longer expand as of ~ 06:55 
UT. Correspondingly, the simulation results show that 
the CME has stopped expanding significantly in a lateral 
direction by this time. (The reader is referred to movie 
CORljointl .mov). 

4.4. Three Dimensional Magnetic Field Topology 

Figures [7] and [8] show time-series plots of the simula- 
tion results matched to the STEREO-B (on-disk) and 
STEREO-A (limb) viewing angles, respectively. (The 
reader is encouraged to view the movies that correspond 
to these figures, SA.mov and SB. mov). 

The inner sphere shows the photosphere with the ra- 
dial magnetic field strength from the magnetogram data. 
The outer sphere (light grey) is at height of 1. li?© (70 
Mm) and represents the altitude at which coronal waves 
are observed. The white contours represent the density- 
enhanced front (same as displayed in Figure [3]). The 
green shade represents an iso-surface of mass density 
with a base ratio (ratio between the current image and 
pre-event image) of 1.1. 

Selected core field lines of the magnetic flux rope are 
drawn in Red and some surrounding field lines of a range 
of sizes are drawn in Blue. Where the core flux rope field 
(Red) reconnects with a surrounding field line (Blue), 
the Blue field line changes to Red, indicating the new 
extended connectivity of the core flux rope field. Recon- 
nections between surrounding magnetic field lines (i.e. 
Blue and Blue) , are shown by the newly reconnected field 
line changing to Yellow. The same magnetic field lines 
have been plotted in both Figures [7] and [5J so that we 
can study the same development from the two different 
perspectives. 

Referring to Figure [3 we see that the green iso-surface 
of increased mass density approximately maps to the 
white contour at each time frame. Figure [5] shows a line 
profile of the density base and running differences, as 
well as the temperature along the path of the coronal 
wave at r = l.li?Q (shown by the black arrow in the 
top panel). It can be seen that the temperature jump 
(indicating the shock front) is ahead of the density in- 
crease associated with the coronal wave. This means that 
the green shade represents the CME front and not the 
shock. Indeed, comparison of the green iso-surface with 
the white-light simulation and observational data (Fig- 
ure [6]) further suggests that it represents the outer-most 
shell of the expanding CME. The existence of a major re- 
connection (discussed below), further suggests that the 
green iso-surface represents the actual CME rather than 
a shock, since there is a close matching between the re- 
connected magnetic field line (Red) and the iso-surface. 
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Reconnection of a given field line first occurs when the 
white contour reaches the vicinity of one of the recon- 
necting field lines. This applies both for reconnections 
between the flux rope and surrounding magnetic field 
(Red), as well as for reconnections between surrounding 
fields (Yellow). We therefore conclude that the reconnec- 
tions are directly driven by the expanding CME. 

Further, we observe a major reconnection between the 
core flux rope field (Red) and the overlying field (Blue) 
at 06:05 - 06:10 UT (marked by narrow white arrows in 
Figures [7] and [5]). This reconnection transfers the connec- 
tivity of the core magnetic field from near the equator to 
a latitude of ~ 60°. On comparison with Figure [51 this 
reconnection explains the development of the secondary 
cavity ( £|4.3[) , and dramatic lateral expansion of the CME 
in the low corona at this time. 

These results show that where the core magnetic flux 
rope is able to reconnect with overlying or surrounding 
favorably orientated field, the CME "opens" up the low 
corona to a very large lateral extent. A clear example was 
report ed in an observational study by iManoharan et all 
( 1996) concerning reconnection with a trans-equatorial 
loop and the subsequent formation of two dimmings in 
quiet Sun regions. Figures [7] and [5] show that quiet Sun 
loops across the whole range of sizes are deflected by 
the CME expansion. Where the orientation is favorable, 
reconnection occurs. Whether CMEs are large-scale by 
nature or become large-scale throug h nurt ure was con- 
sidered by Ivan Driel-Gesztelvi et al.1 (|2008f ). From our 
results, we conclude that CMEs become large-scale by 
nurture, through stretching of the magnetic field and re- 
connection in the l ow corona with t he surrounding mag- 
netic environment ([Pick et al.lll998ft . 

5. DISCUSSION 
5.1. Physical cause of the bright front emission 

It is difficult to separate temperature and density 
effects within single narrow bandpass observations (such 
as 195 A used by EUVI) because line-of-sight effects 
make the optically-thin EUV data complex to interpret, 
and because intensity changes can theoretically be 
caused by alterations in temperature and/or density. 
Indeed, observations have co ntributed evidence in favor 
of b ot h temperature (e. g . IWills-Davev fc T hompson 
19991: IGopalswamv et al.l 120001) and densit y (e!g7 
Warmuth et all 120051: IWhite fc Thompson! [20051 : 
Wills-Davcy 2006) enha ncements. 

IDelannee et al.l (|2008l ) find that the plasma can be 
brightened from both a current shell around the expand- 
ing flux rope via Joule heating, as well as from an in- 
crease in density due to compression. They find that the 
emission from the current shell generates a more conspic- 
uous brightening than that from the plasma compression; 
however, they note that "the dissipation of the current 
densities at low altitude would not be responsible for 
the observed structure" . This work is consistent with 
our simulation results presented in Figure^ which show 
that the coronal wave bright front is well described by 
the mass density enhancement at a height of l.li?©. We 
note that this does not exclude some additional contri- 
bution from temperature effects. Indeed, some heating 
component is expected as the plasma is compressed. 



5.2. What causes the enhanced mass density? 

Delanne e fc Aulanierl (|1999l ) conjectured that fast ex- 
pansion of the magnetic field should compress the plasma 
at the boundaries between expanding stable flux do- 
mains, naturally leadin g to the enhanced emission. 
IChen et all (|2002l . 12005ft showed that stretching of the 
overlying magnetic field leads to compression of plasma 
in the legs of the CME, pro ducing an intensity enhance- 
ment. lAttrill et all (|2007al ) suggested that magnetic re- 
connections between the expanding CME and favorably 
oriented surrounding quiet Sun magnetic field drive weak 
flare- like brightcnings making up the bright front. The 
simulation demonstrates that the plasma is indeed com- 
pressed as a result of the expansion of the CME (Figure 
[3]) and that stretching of the overlying magnetic field 
and magnetic reconnection with surrounding field do oc- 
cur (see Figures [7J and [5]) . Although we only display se- 
lected field lines from the simulation (thus probably miss- 
ing some stretching and magnetic reconnection events), 
these mechanisms are still constraine d to the footpoint s 
of overlying magnetic field lines (see IChen et alj [2002) , 
and locations for favorable reconnection. They are not 
responsible for the emission of the entire bright front, 
only to the higher concentrations of intensity within it 
(Figure El). 

In order to directly drive compression between the ex- 
panding CME and surrounding magnetic field to the spa- 
tial extent covered by the coronal wave, it is necessary 
that the CME really expands to a massive, global extent 
low down (< 200 Mm) in the corona. Our simulation re- 
sults confirm (Figure [6] and ^4.4[) that this is indeed the 
case. 

5.3. Persistent Brightenings 

Stationary brightenings have previou sly been reported 
at coronal hole boundarie s by e.g. iThompson et all 
(fl99ll : lVeronig et al.l (|2006h : lAttrill et al.l (|2007bf ) and at 
separatrices formed in the l arge- s cale magnetic topolog y 
(e.g. IDelannee fc Aulanierl IT9991 : IDelannee et~aTl 120071) . 
Brightening "A" in Figure [1] is likely due to compression 
between the expanding CME and "open" magnetic field 
of the coronal hole to the East of AR 1012, since this 
brightening is also seen in the simulation results (Fig- 
ure [3]) , which show enhanced mass density. Brighten- 
ings "B" and "C" in Figure Q] are located at the same 
place as ongoing magnetic reconnec tions seen in t he sim- 
ulation (Figure [7] and Figure [TO]). lAttrill et all (|2007aft 
reported persistent brightenings similar to these at the 
outer edge of deep, core dimming regions for two events 
during solar minimum of cycle 23. They suggested that 
these brightenings may be due to ongoing reconnection 
between the expanding flux rope and surrounding, favor- 
ably orientated magnetic field. Our results for this event 
are consistent with such an interpretation. 

5.4. Secondary Dimmings 

Behind the expanding bright front, we detect local- 
ized regions of secondary dimming (Figure [4]). Secondary 
dimm ings were originally reported by [Thompso n et al.1 
(2000). Such dimming s may be und erstood in a wave 
context (e.g. as in IWu et all [20011) . since a rarefac- 
tion sh ock develops trailing a l arge-amplitude pertur- 
bation ([Landau fc Lifshita fl987). However, such dim- 
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mings would be short-lived, with a duration on Alfven 
timescales (~ few minutes) contrary to observations 
(|Cliver et al.ll2005HDelannee et al.ll2007t ). Although sec- 
ondary dimmings have a much lower average intensity 
level (e.g. ~ 50 counts/pixel) before the event, compared 
to the core dimming (~ 100 counts/pixel) , the relative 
magnitudes of both the core and secondary dimmings are 
substantial (e.g. ~ 50% and ~ 20%, respectively). Fur- 
ther, like the core dimmings, the secondary dimmings 
remain at a reduced intensity level for an extended pe- 
riod (> 1 hour). 

The locations of these secondary dimmings also ap- 
pear to be closely associated with the magnetic fields 
reconnected through CME expansion. For example, Fig- 
ure [6] shows a secondary dimming to the north of the 
source region that corresponds to the secondary cavity in 
the simulation seen at 06:05 UT. This dimming extends 
the CME cavity northward (see CORljoinl .mov). Fur- 
ther evidence of this can be found in Figure QT] (also see 
iMandrini et al.ll2007f ). These results show that overlying 
and neighboring magnetic field is "opened" through mag- 
netic reconnection, extending the CME footprint in the 
low corona. Therefore, we also interpret the secondary 
dimmings in this event as being due to density deple- 
tion, although spectral diagnostics have yet to confirm 
or refute this interpretation. 

5.5. Understanding the two-component bright front 

Our results show that the bright front observed in base 
difference EUV data and the CME are strongly coupled 
(e.g. Figure [5]). These higher- intensity brightenings are 
due to the CME compressing the plasma (against both 
surrounding and overlying magnetic field). The brightest 
concentrations in the data and simulations in Figure [3] 
show correspondence with the regions of reconnection in 
Figure El both red and yellow field lines. Figure [10] shows 
a direct comparison. 

When the CME has expanded to its maximum lateral 
extent, the brightest parts of the coronal wave either dis- 
appear or become stationary before fading ( 34.2. ip . This 
is the result of multiple factors: the CME is no longer di- 
rectly compressing plasma, the overlying field has already 
stretched, and magnetic reconnections with surrounding 
favorably-orientated field have had time to occur. 

What remains is a weaker, more uniform component 
that is consistent with an MHD wave interpretation. The 
dynamic expansion of the CME is a highly energetic, im- 
pulsive event, therefore wave(s) are expected to be gen- 
erated. The simulation results show that this weaker 
component exists throughout the expansion of the CME. 
The later frames of the simulation show that it contin- 
ues to expand even after the considerable CME lateral 
expansion has finished. In this later stage, the cor onal 
wave is freely propagating (e.g. iVeronig et al.ll200"8l ). In 
the running difference EUVI-A data (195ajrdif f .mov), 
which highlights the leading edge of the disturbance, the 
Western expansion can be followed considerably later 
than in the base difference images, until at least 06:55 
UT (c.f. 06:25 UT, the reader is encouraged to com- 
pare the running and base difference movies for EUVI-A: 
195a_rdif f .mov and 195a_dif f .mov). It is more likely 
that the weaker component can be detected in running 
difference images, which better show subtle changes. 

With this analysis, we believe we are able to recon- 



cile the different (wave and non-wave) interpretations 
of coronal waves. When EIT waves were first discov- 
ered, they were studied using running difference data. 
This method highlights the (often faint) leading edge of 
the disturbance, making it useful for identifying waves; 
this technique was probably perpet uated because re- 
searchers (e.g. I Thompson et al.l 119990 originally identi- 
fied their observations as a strong candidate for the pre- 
dicted ccrojiajxcarnterp^rtto i Hiechromospheric Moreton 
wave (|Moretonlll960l : IUchidalll968h . 

In the late 1990s, Delannee et al., (and later Chen 
et al., and Attrill et al.), preferentially used base dif- 
ference images because they show real bri ghtenings and 
dimmings fe.g. iChertok fc Grechnevl l2005). The motiva- 
tion for using base difference images is due to the focus 
of these authors on coronal dimmings, which are strongly 
connected with coronal waves and CME events. It is not 
possible to study coronal dimmings with running differ- 
ence images. However, base difference images do not 
show faint features so well. We have shown that the base 
difference brightenings are closely linked to the CME and 
magnetic field evolution; hence, the development of non- 
wave models. 

For some time, both of these methods have been ap- 
plied, with different studies producing disparate con- 
clusions. In most cases, however, researchers have at- 
tempted to find a single solution-either wave or non- 
wave-applicable to all aspects of diffuse coronal wave 
events. Over time, this has led to seemingly contra- 
dictory evidence, selectively supporting wave or non- 
wave models, depending of the focus of the study. 
iZhukov k. Aucheri (|2004D first introduced the concept of 
a coupled coronal wave, consisting of an eruptive mode 
and a wave mode, based on comparative analysis of EIT 
running and base difference data. Our results are con- 
sistent with such a picture. The combined observational 
and simulation results presented here allow us to firmly 
establish and understand the contribution from each of 
the various mechanisms. In retrospect, it is not surpris- 
ing that wave and non-wave interpretations have failed 
to be reconciled when the different data sets highlight 
different things! 

5.6. Implications for impulsive electron events 

In interplanetary space, field lines following the Parker 
spiral connect the western longitudes of the Sun to space- 
craft at 1 AU. When an impulsive electron event occurs 
in the corona, energetic particles travel along these field 
lines to 1 AU. However, sometimes impulsive electron 
events are clearly related to flares that occur on the 
eastern half of the sola r disk, up to ~ lJ?p from the 
Parker spiral footpoint. Kruckcr ct al. (1999) suggested 
that EIT waves might explain how the flare site and Sun- 
spacecraft magnetic field lines are connected. They con- 
cluded that at the time of electron release, the EIT wave 
had not expanded far enough to reach the footpoints of 
the Parker spiral. However, by considering the EIT wave 
as a fast-mode MHD wave (which moves faster at higher 
altitudes due to decreasing density), they showed that 
the calculated wave front at higher altitude (~ 1.5Rq) 
is fast enough to connect the flare site with the Sun- 
spacecraft field line. 

We would like to suggest an alternative possibility. Our 
simulation results (Figures [7] and [8]) show that recon- 
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nection between the core magnetic flux rope (intimately 
connected with the flare region) and the overlying or sur- 
rounding magnetic field occurs at a similar height range. 
The reconnection(s) subsequently displace the flux rope 
connectivity out of the flare region to a distance of ~ 
IRq, consistent with observations. 



6. CONCLUSIONS 

Using a combination of multi- wavelength observational 
analysis and a global MHD simulation driven by real 
data, we study the 13 February 2009 coronal wave - CME 
- dimmings event observed by STEREO, in quadrature. 
We find that the diffuse bright front emission is primar- 
ily due to mass density enhancement. This is caused by 
a combination of both wave and non-wave mechanisms, 
both of which are directly-driven by the CME, which ex- 
pands to a considerable lateral extent in the low (< 200 
Mm) corona. The reorganization of the magnetic field 
through reconnection facilitates lateral expansion of the 
CME, leading to far-reaching compression and "open- 
ing" of the surrounding magnetic field. This is further 
evidenced by secondary dimmings that form across the 
solar disk. 

We find the diffuse coronal wave front displays a dual 
structure, consisting of a brighter and a weaker compo- 
nent. The brighter component, observed primarily in 
base difference EUVI data, is due to plasma being com- 
pressed by the expanding CME (against both surround- 
ing and overlying magnetic field). Some of the bright 
patches correspond to regions of reconnection where the 
field orientation is favourable for this to occur. This non- 
wave component maps directly to the CME footprint at 
every stage of the evolution. Thus, there is a strong 
coupling between the development of the coronal "wave" 
bright front, CME and associated dimmings. 

The weaker component, which is most likely to be 
detected in running difference EUVI data is present 
throughout the event, and ultimately decouples from the 
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bright component after the CME ceases lateral expansion 
late in the event. This suggests that the weaker compo- 
nent is an MHD wave, initially driven by the expanding 
CME, later becoming freely propagating. 

This work demonstrates the considerable insight 
gained from advanced numerical simulations well con- 
strained by observations. We hope that this work can 
progress the coronal wave community away from divisive, 
separatist theories toward a more cohesive, holistic ap- 
proach to understanding the complexity of EUV coronal 
waves. Future work should focus on the combined anal- 
ysis of other well-observed events and on what coronal 
"waves" can tell us both about their driving CMEs, and 
the structure and dynamics of the surrounding corona. 
The potential for coronal seismology can now be pur- 
sued with the confident identification of the true wave 
component. 

As this study demonstrates, detailed global MHD sim- 
ulations are essential for furthering development of com- 
prehensive physical models. We must now include rigor- 
ous quantitative data analysis for comparison with such 
models. This goal will be forwarded both by the contin- 
uing development of automated measurement techniques 
and by the upcoming launch of the Atmospheric Imaging 
Assembly aboard the Solar Dynamics Observatory. 
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Fig. 1. — Top panels show pre-event EUVI-B 195 A Fe XII images at 05:15 UT. Black contours are overlaid, outlining regions of low 
intensity. Bottom panels show base difference EUVI-B images at 06:05 UT. The same contours as in the top right panel are overlaid on the 
base difference data. The persistent brightcnings (white regions) lie at the boundaries of the deep, core dimming regions. A low-latitude 
coronal hole lies to the east of brightcnings "A" and "B" . Brightenings marked "A" are evident in the simulation results as well (see Figure 
[3t . and hence are due primarily to a density increase. Brightenings "B", and "C" on the other hand, are much more closely linked to 
regions of ongoing magnetic reconnection between the core flux rope magnetic field and surrounding, favorably orientated magnetic field 
(see Figure (7JI. 
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Fig. 2. — Comparison of Stereo-A COR1 images and synthetic white-light images from the simulation at 5:55, 6:15, and 6:35 UT. Each 
row shows a COR1 base difference image (left), base difference image of the simulated white-light scattering (center) and a COR1 running 
difference image (right). The white arrows in the right panels indicate the maximum lateral extent of the CME as seen in the COR1 running 
difference movie - the reader is encouraged to view the COR1 running difference movie accompanying this Figure (C0Rl_rdif f .mov). The 
expansion of the subtle outermost CME shell to the North is more easily discerned than the expansion to the South, where streamers 
complicate the observations. See also the corresponding movies: CORljointl .mov and C0Rl_Joint_Zoom.mov. 




EUVl(B) 1 95 Q6:05 UT 







Fig. 3. — A time series for the period 5:25-6:35 UT comparing STEREO-B EUVI Fe XII observations and the simulation results. Each 
pair compare STEREO-B EUVI base differences image (left) with base differences of simulated mass density at height of l.lRs u „ (right). 
Overlaid on the base difference images arc dashed white circles which act as a guide and are drawn by eye, to indicate the maximum extent 
of the bright front. The movies for this figure are: 195b_dif f .mov and densityf rontdif .mov. The corresponding data from the viewpoint 
of STEREO-A arc included as movies: 195a_dif f .mov and densitysidedif .mov. 
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Fig. 4. — Output from the automatic dimmings extraction algorithm (Attrill & Wills-Davey, 2009). The deep, core dimming is extracted 
in the immediate vicinity of the active region (c.f. bottom panels, Figure[T]l. In addition, widespread, secondary dimmings are also detected, 
spread across a large fraction of the solar disk (c.f. Figure [3) l. 
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Fig. 5. — The steady-state values of number density (top left), magnetic field strength (top right), temperature (middle left), plasma j3 
(middle right), sound speed (bottom left), and Alfven speed (bottom right) at height of r = I.IRq. 
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Fig. 6. — Left panel shows the COR1-A running difference image at (06:05 - 05:55) UT. Center panel shows the simulated white-light base 
difference image at (06:05 - 05:35) UT. Right panel shows the EUV1-A running difference image at (06:05 - 05:55) UT. All images are scaled 
to the same size. This figure demonstrates the importance of the role of the simulation in developing an understanding of th e true lateral 
exten t of the CME in the low corona. Fitting the CME as observed in the COR1 data alone (e.g. see Figure 4, [Patsourakos & V ourlidasl 
2009), a large part of the CME in the low corona is missed. When proper consideration of the CME extent in the low corona is made, the 
lateral extent of the CME maps to the coronal wave observed in EUV1 data. The white arrow in the center panel indicates a region which 
is also part of the CME cavity (c.f. bottom panels, Figure[2]at 06:35 UT). 
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Fig. 7. — A time series (from 05:25 - 06:20 UT, for frames until 07:15 UT see movie SB.mov) of the simulation results matched to 
STEREO-B viewing angle. Inner sphere shows the photosphere with the radial magnetic field strength (magnetogram data). Outer sphere 
is at a height of I.IRq colored with white line contours of mass density base difference (same as Figure|3jl. The light Green shade represents 
an iso-surface of mass density with a base ratio of 1.1. The flux rope core field lines are drawn in Red and surrounding field lines are drawn 
in Blue. A Blue field line is changed to Red if it reconnects with a flux rope core field line, and to Yellow if it reconnects with a surrounding 
field line. This figure is also included as a movie: SB.mov. 




Fig. 8. — Same as Figure [7] but matched to STEREO-A viewing angle. (See movie SA.mov). 
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Fig. 9. — An extraction of density base and running differences (second panel), and temperature (bottom panel), at r = l.l-R© along the 
line shown in the top panel. Red line marks the location of the coronal wave front, while yellow line marks the location of the CME shock. 
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Fig. 10. — Top left panel shows a base difference image of mass density at 1.1 R© from the simulation. Top right panel shows the EUVI-B 
base difference data at 05:55 UT. Bottom panel shows an overlay of the mass density and magnetic field from the simulation. Orange 
arrows indicate the different mechanisms that contribute to the bright front observed in the base difference data at 05:55 UT. 
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Fig. 11. — Top left panel shows a base difference image of mass density at 1.1 R0 from the simulation. Top right panel shows the core 
and secondary dimmings at 07:05 UT, extracted using the automatic dimmings algorithm. Bottom panel is an overlay of the extracted 
dimmings and selected magnetic field lines showing the correspondence between reconnected magnetic field (red and yellow lines) with the 
locations of secondary dimmings at 07:05 UT. 



